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Abstract 



We consider a gravity dual description of time dependent, strongly interacting large--/V c 
M = 4 SYM. We regard the gauge theory system as a fluid with shear viscosity. Our 
fluid is expanding in one direction following the Bjorken's picture that is relevant to 
RHIC experiments. We obtain the dual geometry at the late time that is consistent 
^ ■ with dissipative hydrodynamics. We show that the integration constants that cannot be 

determined by hydrodynamics are given by looking at the horizon of the dual geometry. 
Relationship between time dependence of the energy density and bulk singularity is also 
discussed. 



*E-mail: nakamura@ihanyang.ac.kr 
^E-mail: sjsin@hanyang.ac.kr 



1 Introduction 



One of the attractive aspects of AdS/CFT k is applicability to the real systems after certain amount 
of deformations. In fact it has been suggested that the fireball in Relativistic Heavy Ion Collision 
(RHIC) can be explained from dual gravity point of view [31313101) since the quark-gluon systems 
created there are in the strong coupling region [2] . Although the YM theory described by the standard 
AdS/CFT is large- N c M = 4 SYM theory, there are many attempts to construct models closer to 
QCD 7 . SUSY is not very relevant in the finite temperature context since it is broken completely. 

Since the RHIC fireball is expanding, we need to understand AdS/CFT in the time dependent 
situation. Recently, Janik and Peschanski [HI [H] discussed this problem in non-viscous case. They 
use the conservation law and conformal invariance together with the holographic renormalization 
|10| II lj to express the bulk geometry with given boundary data. As a result, the bulk geometry 
reproduces the basic features of Bjorken theory It is a ls° pointed out that inclusion of shear 
viscosity, although the value is small, is very important in the analyses of real RHIC physics since it 
plays an essential role in the elliptic flow (see for example, ^3 El)- In fact, the shear viscosity at 
the strong coupling limit was calculated for the M = 4 SYM systems in Ref. ^5] using AdS/CFT. So 
it is natural to ask how the bulk geometry change if we include the viscous effects in the boundary 
theory. 

In this paper, we establish the dual geometry in the presence of shear viscosity by using the hydro- 
dynamics as the boundary data. Although our gauge theory is not QCD, we hope there is universal 
features in the character of strongly interacting gauge theory systems. In fact hydrodynamics, which 
is our input does not ask much about the details of the microscopic particles and the interactions 
once the equation of state is given. Therefore we have a chance to extract useful information on the 
macroscopic properties of the real quark-gluon fluid based on this universality. 

If what we get is consistency with fluid dynamics, there would not be much point to consider 
AdS/CFT dual of it. In fact, the holographic dual of the hydrodynamics contains much more informa- 
tion than the hydrodynamics since AdS/CFT already contains essential information of microscopic 
gauge theory dynamics. For example, we will show that the holographic dual of the hydrodynamics 
gives integration constants in the hydrodynamic equations that cannot be determined by hydrody- 
namics alone. The dual geometry also gives a simple derivation of Stefan-Boltzmann's law in strongly 
coupled regime with precise Stefan-Boltzmann constant. 

The organization of the paper is the following. In Section [21 we analyze time dependence of 
the system in the framework of the relativistic hydrodynamics. We also review the basics of the 
dissipative relativistic hydrodynamics in order to clarify our setup. Section [21 gives the analysis in 
the gravity dual. We review the basic framework of the gravity dual and present some results for 
non- viscous cases obtained in Ref. [S]. The main results of the present work will be given in Section 
IH. HI where the late time dual geometry is proposed and consistency with the hydrodynamic analyses 
is checked. We will show that the holographic dual of the hydrodynamics contains more information 
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than the hydrodynamics. We also make comments on the regularity of the bulk geometry in Section 
@J We conclude in the final section. 

2 Relativistic hydrodynamics with shear viscosity 

We begin with a short review of the relativistic hydrodynamics with dissipation. For the relevance 
to the RHIC fireball, we assume that it is described by the finite temperature theory of a variant of 
N = 4 SYM. We also follow the Bjorken's picture In the Bjorken's model, the system undergoes 
one-dimensional expansion (Bjorken expansion) along the collision axis of the heavy ions, and the 
fluid of the quarks and gluons has boost symmetry in the so-called central rapidity region |12j . We 
shall consider only the late-time regime of the Bjorken expansion where the time evolution is slow 
enough to employ approximations. 1 

The energy-momentum tensor in the framework of relativistic hydrodynamics is known to be 2 

T^ u = {p + P)u^u v + Pg^ +t» u , (2.1) 

where p, P are the energy density and the pressure of the fluid, and = (7, jv) is the four- velocity 
field in terms of the local fluid velocity v. t^ v is the dissipative term. In a frame where the energy 
three-flux vanishes, r^ v is given in terms of the bulk viscosity £ and the shear viscosity r/ by 

= - v (A^ x V x u u + A^Vau" - \ V x u x ) - £ V x u\ (2.2) 

under the assumption that t^ v is of first order in gradients. We have defined the three- frame projector 
as A"" = g^ + vPu". 

In this paper, we consider pure M = 4 SYM theory whose energy-momentum tensor is traceless. 
Now, the trace of the energy-momentum tensor is given by 

r; = - P + 3P-3ev AU A . (2.3) 

Demanding Tff = for all the possible frames where (|2.2|) is valid, we obtain 

£ = 0, and p = 3P. (2.4) 

Notice that the bulk viscosity in the realistic RHIC setup might also be negligible. (See for example, 
Ref. Q3|.) 

We assume that our fluid system is boost invariant following Bjorken since it is actually 

supported by experiments. We want to take a "co-moving frame" where each point of the fluid labels 

the coordinate, a concept called Lagrangian frame in fluid dynamics. In this frame all the fluid points 

1 Realistic model should contain three-dimensional expansion as tried in Ref. [5j. There, it was suggested to use a 
dual of three-dimensional cosmic expansion. 

2 The convention of the signature of the metric is ( — , +, +, +) in this paper. 
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are at rest by definition, hence all the fluid points share the same proper time. We can use the rapidity 
of each fluid-point as a spatial coordinate and the common proper time of each fluid-point as a time 
coordinate. Therefore a local rest frame (LRF) of the fluid can be given by proper time(r)-rapidity(y), 
whose relationship with the cartesian coordinate is (x ■ X • X • X 3 ) = (rcoshy,rsinhy,2; 2 ,x 3 ). We 
have chosen the collision axis to be in the x 1 direction. 

The Minkowski metric in this coordinate has the form of 



ds z 



-dr 2 + T 2 dy 2 + dx 2 



(2.5) 



where dx\ = (dx 2 ) 2 + (dx 3 ) 2 . We assume that the collision happened at r = and we consider 
only r > region. Note that \y\ ~ oo corresponds to the fronts of the expanding fluid. Therefore, 
the whole region on the y-coordinate axis is occupied by the fluid. We also assume that the fluid is 
extended in the x 2 , x 3 directions homogeneously. Since the real fireball produced by RHIC experiment 
is localized, the set up we use is an idealized one. Nevertheless, the present setup is proper since we 
are interested in the central rapidity region. 



The four-velocity of the fluid at any point in the LRF is 
energy-momentum tensor to be diagonal: 



(1,0,0,0), and this makes the 



I p 


V o o 



111 

3 T 






p + 11 









P + In 

r + 3r 



(2.6) 



We have three independent quantities, p, P and r\ in (|2.6|) . However, the energy-momentum conserva- 
tion, V^T^ = 0, together with the equation of state p = 3P, reduces the number of the independent 
quantities to be one. One finds that the energy-momentum tensor is written by using only p in the 
following way: 



rpflU 



( P 








7? irP-rp) 








p+\rp 










\ 



(2.7) 



P + 2 T P J 



where p = By identifying (|2.6|) with ()2.7|) . we obtain the following differential equation that 
connects 77 and p: 3 



dp 4 p 4 i| 
dr 3 r 3 r 



+ dr ( 2 - 8 ) 



3 The the equation 12. 81 turns out to be the same as the one appearing in so-called first order (or standard) dissipative 
relativistic hydrodynamics. (See for example, Ref. |16j and the cited therein.) It is known that the first order formalism 
has a problem of acausal signal propagation. However, it gives good enough results for our purposes. The details of 
the causal dissipative relativistic hydrodynamics and consistency of our analysis are shown in Appendix 151 
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Note that both of p and r] depend on the proper time r in general. Let's assume that the shear 
viscosity evolves by 

, = (2.9) 
where 770 is a positive constant. The solution of (|2.8j) is then given by 

p(T) = ^ + r=h^ (for " * 1/3) ' (2 - 10) 

where is a positive constant. For /3 < 1/3 case, the viscous corrections in the hydrodynamic 
quantities become dominant in the late time, which invalidates the hydrodynamic description. If 
P > 3, the shear viscosity term is sub- leading in the late time behavior as we expect. Therefore we 
will consider only j3 > | case from now on. 

The proper time dependence of the temperature T can be read off by assuming the Stefan- 
Boltzmann's law p oc T 4 : 

T = T °{^lTh-^ + -)- < 2 ' 12 > 

In the static finite temperature system of strongly coupled N = 4 SYM theory, it is known that 
77 oc T 3 ^3]. Let us assume that the same is true in the slowly varying non-static cases. Then we set 
(3 = 1: 

r,= ™. (2.13) 
r 

We know p ~ T 4 and 77 ~ T 3 cannot be consistent without an additional term in (|2.13j) . but the 
correction term is negligible in our case. One can check the consistency of our approach in Appendix 
IO The temperature behavior is then given by 

T = r »(^-^7 + -)- < 2 - 14 > 

We can evaluate the entropy change in the presence of shear viscosity by using hydrodynamics. 
The conservation of energy-momentum tensor can be rewritten as 

¥ + *=£ < 215 » 

Using the nature of the one-dimensional expansion, the above can be rephrased as 

T d(rs) _ 4r? 
dr 3r ' 

where s denotes the entropy density and rs = S is the entropy per unit rapidity and unit transverse 

area. 4 Notice that in the absence of viscosity, S is constant. Now, the entropy per unit rapidity and 

4 A precise definition of S is the entropy within a unit 3d region on the (y,x 2 ,x 3 ) coordinate. The voiume of this 
region is r and it is expanding with time in the x direction. 



(2.16) 
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unit transverse area has time dependence, 

4// 



S{t) = J dr 



3tT 

Z?° -2/3 , n u-i/3^ 



= Soo -2^t-^ + 0(t-^), (2.17) 

due to creation of the entropy by dissipation. However, the creation rate of the entropy slows down 
with time. Soo, that is the entropy per unit rapidity and unit transverse area at r = oo, is an 
integration constant which we cannot determine in the framework of hydrodynamics. We will show 
that its precise value is given by using the gravity dual in Section [3 



3 Holographic Dual of Hydrodynamics 

In this section, we will find a five-dimensional metric which is dual to the hydrodynamic description 
of the YM fluid in the previous section. The basic strategy is to use the Eistein's equation together 
with the boundary condition given by the energy-momentum tensor at the boundary |10Ull| l8j. We 
consider general asymptotically AdS metrics in the Feffer man- Graham coordinate: 

2 2 9^dx^dx v + dz 2 
ds=r^ = , (3.1) 

where x^ = (r, y, x 2 ,x 3 ) in our case. r$ = (<i'Kg s N c a' 2 ) l l 4: is the length scale given by the string 
coupling g s and the number of the colors N c . The four-dimensional metric g^ v is expanded with 
respect to z in the following form |10| 111): 

9flv (r, z) = 0<g>(r) + z 2 g$(r) + z 4 g$(r) + z 6 g$(r) + ■■■. (3.2) 

gffl is the physical four-dimensional metric for the gauge theory on the boundary, that is given by 
(|2.5[) in the present case. The g^L ■ 's depend only on r because of the translational symmetry in the 
x 2 , x 3 directions and the boost symmetry in the y direction in our setup, gffl is found to be zero. 
We can identify the first non-trivial data in (|3.2|) . g$ , with the energy-momentum tensor at the 
boundary 10 : 

$ = ^<^>, (3.3) 
r o 

where G5 is the 5d Newton's constant given by G5 = 87r 3 a' 4 g 2 j 'r\ in our notation. For the time 
being, we set ^ttG^ = 1 and ro = 1. The higher-order terms in (|3.2j) are determined by solving the 
Einstein's equation with negative cosmological constant A = —6 |1()| IHj: 

1 

Rmn — -^GmnR — 6Gmn = 0, (3.4) 

where the metric and the curvature tensor are for the five-dimensional ones of (|3.1|) . g^J 1 ^ is described 
by gffi 2 ^ , g^v , ■ ■ ■ , gjS) through solving the Einstein's equation. In other words, we can obtain 
the higher-order terms in ()3,2j) recursively by starting with the initial data gffi (~Minkowski) and 
9$ (~ 
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3.1 Static cases 



In order to demonstrate the use of the above procedure, we first workout the static case using cartesian 
coordinate = [t : X , X . X i ) instead of the Bjorken's coordinate. The energy-momentum tensor in 
this case is given by 



= diag 0, p/3, p/3, p/3) 



(3.5) 



The result of the above procedure gives the solution of the Einstein's equation in Fefferman-Graham 
co-ordinate: 



ds 2 = 1 1 J 1 i z4 } 2 dt 2 +(i+ |/) (dxj + dxi + dxi)\ + d -i. 



z" 



l+§^ 



(3.6) 



which is equivalent to the AdS-Schwarzschild Black hole. 5 The Hawking temperature is given by 



T H = V2/(z tt 



(3.7) 



where zq(t) = [3//o] 1//4 is the position of the horizon. By restoring AttGq and ro = (4:7rg s iVca' 2 ) 1 / 4 , 
and by identifying Tjj with gauge theory temperature T, we obtain the Stefan-Boltzmann's law 6 



(3.8) 



This result agrees with Ref. [T7| . 

3.2 Non-viscous time dependent cases 

Coming back to Bjorken's setup, Janik and Peschanski obtained the late time bulk metric [S] by the 
above procedure. Here we briefly review their work in our language. By using (|2,7|) and ()2.1U|) . the 
energy-momentum tensor for non-viscous case is explicitly written as 



( 


PO 

7373 











\ 







T 3r 4 /3 


















PO 
3r 4 /3 







V 











PO 
3r 4 /3 


J 



(3.9) 



Then the metric is given by 



-1 + — —z +<J(z ), 



r 4/3' 



9yy 
r 2 

9xx 



1 + 
1 + 



3r 4 / 3 
Po 



z A + 0{z & ), 
z A + 0(z 6 ), 



(3.10) 



3r 4 / 3 ' 

5 Notice that the met ric is map ped to the standard form of the AdS-Schwarzschild metric through the coordinate 
transformation z = zj ' \/l + z 4 /zq. 

6 Actually Stefan Boltzmann's law is about intensity I = cp/A in terms of temperature. But we use the terminology 
abusively to name (13.81 . 
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where g xx = 522 = 533- Notice that our Minkowski metric is given by (|2.5|) , Let's focus on the 
late time behaviour of the metric, since T^y is given only for the late time.''' If we take t — ► 00 
limit naively, what we obtain is just the Minkowski metric (|2.5|) . To extract a non-trivial result, the 
authors of Ref. [B] take the limit such that g^> z 4 does not go to zero nor infinity as r — > 00: 

t — ► 00 with — = v fixed. (3-11) 
j- 1/6 

Then by solving the Einstein's equation recursively up to certain order of z, g TT , g yy /T 2 and g xx have 
the following structure: 

f {1) (v) + f i2) (v)/r^ + .... (3.12) 
By neglecting the 0(r -4 / 3 ) quantities, they obtained an analytic expression of the late time metric: 

d - 2 = ?{- ".; y^'' * 3 + f 1 + t^j) (tW + dii) ) + f • (3 - i3) 

We can explicitly check that the above metric satisfies 

G lm {Rmn - \g MN R - QGmn) ~ 0(l/r 2 ). (3.14) 

Notice that (|3.13() is a black hole in AdS space with time-dependent horizon. The time dependence 
of the entropy and the Hawking temperature from the metric (|3.13[) reproduces the Bjorken's results 
|12j S ~ constant as well as T ~ t -1 / 3 . Interestingly, they observed that the regularity of the 
geometry at the horizon uniquely select the power of time evolution of energy density p ~ r -4 / 3 , 
which is a consequence of hydrodynamics. (See Section 0] for the details.) 

If we replace p = ^ with p = ^ ; we f] nc j that the left-hand-side of ()3.14j) is at the order 
of 1/(t 2 logr). In this sense, this replacement makes another late time solution of (|3.4|) . We will use 
this solution in Section |1J 



3.3 Viscous cases 



Let us come back to our main interest to obtain the bulk geometry in the presence of shear viscosity. 
The energy-momentum tensor for = 1 is written by using 1)2.7(1 and (|2.1Uj) as 



/ Pa 

I T 4/3 



- )1U 







V 









( pp. _ 230 A 
V37473 T^J 









PQ, 




\ 






3r 4 /3 



(3.15) 



/ 



The metric components, g TT , g yy /T 2 , g xx have the following structure by solving the Einstein's 
equation recursively: 



fW(v) +r]0 h^(v)/T^ + p\v)/r^ + ..., 



(3.16) 



The slow time evolution is necessary to justify the hydrodynamic treatment of the fluid. 
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Note that the viscosity dependent terms exist at the order of r~ 2//3 and these are more important 
than the higher-order terms neglected in \3.13\) . We are considering the late time region t>1. But 
to see the effects of viscosity, we need to keep the terms at least to the order of t~ 2 / 3 . In this paper 
we consider the viscosity effects to the minimal order. 

Now we solve the Einstein's equation recursively. The power series that appear in the solution 
can be re-summed to give a compact form of the metric. For the detail, see Appendix [3 The late 
time 5d bulk geometry is given by 



1 



z 2 




-2t 



ds Z = ^<~- ^^+11 + ^)1^—^1 r 2 d 2 y + 




dz 2 

+%, (3-17) 

z^ 



where 



and p = -rpr -. (3.18) 



' ~ p T 2 / 3 y r 4 / 3 r 2 

Notice that the energy momentum tensor (|3.15j) can NOT be written in terms of the whole p{r). It 
is truly amazing that the final metric nevertheless can be written in terms of p(r) (apart from the 
power) in the compact way. This implies that the position of horizon can be determined solely by 
the energy density. 8 

The Hawking temperature in the adiabatic approximation is given by T(r) = \^2/t:(zq(t)), where 
z o( T ) = [3//o(t)] 1 / 4 is the time dependent position of the horizon. Just as the static case, we obtain 

p = ^n 2 N 2 T\r), (3.19) 

by restoring ^ttG^ and ro- The entropy per unit rapidity and unit transverse area is given by 

1 2V2rr 3 
4G 5 zl{r) 

One remarkable thing is that the value of S at r = oo, that cannot be determined by hydrodynamics 
alone, is precisely determined to be 

1/4 



in terms of the initial condition pQ. 

8 One should keep in mind that we are looking for the late time geometry; the metric (13.1711 is correct only to the 
order of 7 and the Oipj 2 ) contributions are not unambiguously determined. The representation of 13.171 is chosen since 
it makes the volume of the horizon finite. 
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Let us check consistency of (j3.2U|) and 1)2.17(1 . The normalized entropy-creation rate is given by 

+ 0(t- 7/3 ) (3.22) 



1 dS VO ,-„_-7/3\ 



from the gravity dual and 



S dr Pqt 5 / 3 



Sdr 3 Tq^ooTS/S 
from (|2.17|) of the hydrodynamics. Comparing (|3.22|) and (|3.23|) . we obtain 

4 po 4,pr 



+ 0{t-'I 6 ) (3.23) 



g 

00 3 T 3 T 



(3.24) 



This is nothing but the relationship among the entropy, the energy (per unit rapidity and unit 
transverse area) and the temperature obtained by thermodynamics at r = oo where the system 
reaches thermal equilibrium. 

Before closing this section, let us give a technical remark to clarify the meaning of the late time 
limit (|3.11|) . The readers might want to skip this paragraph at first reading. One finds that the 
higher order terms we have neglected in (|3.12|) and (|3.16|) contain the terms proportional to i; 6 t -4//3 
for example. On the other hand, the leading order terms in g TT contains arbitrary higher power of v. 
Therefore, neglect of the 0(t -4 / 3 ) terms is justified only when v is larger than 0(1). In fact, we can 
justify the limit (|3.11|) when we extract the thermodynamic quantities of the system. Such quantities 
are associated with the horizon of the black hole and the value of v at the horizon is indeed 0(1) 
constant at the late time. Now, we understand why the naive r — > oo limit on the z-coordinate (the 
limit with fixing z to be constant) is not good for our purpose. The position of the horizon grows 
with time: zq ~ r 1 / 3 . If we treat z to be a constant in the r —* oo limit, the region we describe 
becomes infinitely far from the horizon. In other words, neglect of the terms of (z/r 1 / 3 )™ (n > 0) is 
not justified around the horizon. This means, a suitable coordinate is the ^-coordinate rather than 
the z-coordinate to describe near the horizon at r — > oo. 



4 Conditions on energy density and bulk singularity 

In Ref. jH], singularity analysis was used to select a physical metric. Namely, starting with energy 
density (without viscosity) 9 

P=% (4.D 

it is found that the late time bulk geometry is singular except for a special value of /. More precisely, 
(Rmnkl) 2 & t the order of (r)° has singularity at the horizon except for I = |, the value for the perfect 
fluid. In fluid dynamics, this value of I is determined by the conservation law and the equation of 
state. However the bulk metric knows the correct form of the energy density independently |Sj. 
9 The value of I is restricted to be < I < 4 by the positive energy condition for l|2.7|l 8 . 
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Therefore it is interesting to see whether requiring the regularity of the metric at the horizon gives 
further control over the behavior of the viscous term as well. The late time bulk geometry with 
generic value of (3 can be obtained through similar calculations starting with the energy-momentum 
tensor (}2~Tj) with p given in l|2~TUj) or (f2~TTj) . 

If f3 < 1/3, the viscous correction in p is dominant at the late time and p ~ 
This leads to the singular geometry since the proper time dependence of the dominant term in p, 
1/t 1+ P, is not 1/r 4 / 3 . 

If [3 > 1/3, the viscous correction is sub-leading and we should consider {Rmnkl) 2 to the sub- 
leading order. The late time geometry is the same form of (j3,17|) except the following replacement 

T-t'^Apot^ 3 ), (4.2) 

and with p given by (j2. 1U|) . We find that [Rmnkl) 2 has the following structure: 

,16 _i_ on,, ,12 _|_ yjA„,$ _i_ on„, A 
(1 + w 



2 8(5u; lfa + 20^ + 174^ + 20 W 4 + 5) n( ^ /3) „\ 
[Rmnkl) = M , _ 4yt V U(t ), (4.3) 



where 



3/ rV3 V 3 / V 1 - 3/3 

In the absence of viscosity, the above result is reduced to that of Ref . [S] . The first term in the right- 
hand side of (|4.3|) . that contains the viscous sub-leading corrections, is finite. So the consideration 
of singularity does not give any further restriction to the viscosity term. If j3 > 5/3 the corrections 
due to the viscosity-dependence give smaller effects in the metric than the non-viscous 0(r _4//3 ) 
corrections which are already discarded in the late time geometry. So up to our approximation, the 
viscous effect is not visible in this case. 

For (3 = 1/3 case, the viscous corrections are leading order. We find that the metric (|3.13jl . where 
^7373 is replaced with the second term of p in ()2.11j) . gives the late time geometry 10 as mentioned 
in Section I5~2l (Rmnkl) 2 at the leading order is given in the same form as that of ()4.3|) where w 
is w = ^(p/3) 1 / 4 with p given by the second term of (|2.11|) . There is no divergence at the leading 
order, although the value of (3 is not consistent with the hydrodynamics. 

Altogether, our geometry is regular and the viscosity effect is meaningful in the region of 

1/3 < /? < 5/3. (4.5) 

Indeed, = 1 is within this region. This means that in the late time geometry (or the late time 
fireball dynamics), the viscous term gives visible contribution to the dynamics of the fireball. 

10 In this case, the late time limit should be taken by fixing v = z(log r) 1 ^ 4 /r 1 ^ 3 . 



10 



5 Conclusions 



We considered the gravity dual of large- iV c M = 4 SYM fluid undergoing one dimensional expansion 
with account of shear viscosity. We obtained the late time bulk geometry to the minimal order of 
the viscous corrections in the analytic form. We found that our viscous corrections do not break the 
regularity within our approximation. We also found that the time evolution of the thermodynamic 
quantities given by the late time geometry is consistent with the hydrodynamic analyses. 

We saw that the holographic dual of the hydrodynamics contains much more information than the 
hydrodynamics, since AdS/CFT already contains essential information of microscopic gauge theory 
dynamics. For example, the holographic dual of the hydrodynamics gave a simple derivation of 
Stefan-Boltzmann's law in the strongly coupled region with precise Stefan-Boltzmann constant. The 
integration constant in the hydrodynamic equation was also given by looking at the horizon of the 
dual geometry. 

We believe that by probing the resulting geometry, one can extract many of information of strongly 
interacting system in principle. It is important to compute various physical quantities based on the 
obtained geometry ^Jj. We can also consider various extensions of the present work. Inclusion of the 
higher-order corrections, generalization to the systems with chemical potential, consideration of the 
systems with three-dimensional expansion are possible directions. One can also consider the effects 
of the bulk viscosity whose presence violates the conformal invariance. In the real QCD, conformal 
invariance must be broken and including it might be relevant for more realistic account of RHIC 
fireball. We hope that the present work will shed light on AdS/CFT for non-static non-equilibrium 
systems and holographic description of RHIC physics. 
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A Derivation of the metric ( 13.171 ) 

First, we find the following expressions by solving the Einstein's equation recursively to the order of 



v 



,28. 



9tt{t,v) 



1 + 3a - 4a 2 + 4a 3 - 4a 4 + 4a 5 — 4a 6 + 4a 7 H 




(l + (- 2 + 3a - 4a 2 + 5a 3 - 6a 4 + 7a 5 + • • • 



)) + 0(r" 4 / 3 ) 
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= 1+ °-^( 1 + l( 1 + ! + T + T + T + y + -j) + ^ 4/3 )- 

where a = pqv 4 /3. 

The power series in the right-hand sides are re-summed to give the analytic form of the metric: 



4 . , „„„4 , 



9tt{T, V) 



2 

3 



T 2 i+ 3 3r 2/3 I i+ pc^l g I 1 PQV* 1 1 

*.<t,.) = l + «^-^|l-il^l„ g Jl±^]| +0 ( T -n (A.2) 

We can check explicitly that (|A.2|> is indeed the solution of the Einstein's equation that is accurate 
to the order of r -2 / 3 . The important feature of (jA.2|) is that all the terms at the order of t~ 2 / 3 are 
proportional to tjq. ()A,2|) is rewritten in the following forms up to 0(7 2 ) terms: 



/, _ p (l-27)u 4 n2 

te(T ,„) . f 1+ Mj^^ 1 + 7lo Jl±^|| +0(T - n ( , 3) 



3 



or more simply, 



1 + 3 



T, U 



„4 



r 2 ^ 3 

,,4 



1 + 


Po(l- 


-2 7 )t> 4 




3 


1 - 


Po(l- 


-2 7 > 4 




3 


1 + 


Po(l- 


-2 7 > 4 




3 


1 - 


Po(l- 


-2 7 > 4 




3 



-2 7 

+ 0(r" 4 / 3 ), 



4 \ 7 



(A.4) 



The differences among (|A~2l) . (lA~3l) and (TOl) are at the order of 0(t" 4 / 3 ) (that is the same order of 
0(7 2 )) and they share the same terms to the order of 7. 
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B Second order formalism of dissipative relativistic hydrodynamics 



In this appendix, we briefly introduce the causal dissipative relativistic hydrodynamics that is also 
referred to the second order dissipative relativistic hydrodynamics. We shall show that the first order 
formalism we have employed gives a good approximation of the second order formalism at the late 
time in our setup. 

It is known that the first order formalism of the dissipative relativistic hydrodynamics has a 
problem; the viscous and the thermal signal propagates instantaneously and causality is broken. (See 
for example, JH].) I n the second order formalism j!8j . the relaxation time of the fluid is introduced 
to maintain the causality. The energy density evolution equation in our setting is 

? (B.« 

CLT 3 T T 

where 

4 r\ 

= — — (the first order formalism), (B.2) 
d& 4 r] 

T Tr—j- = + (the second order formalism), (B.3) 
and 7V is the relaxation time of the system. Note that the r n — * limit gives the first order formalism. 

Let us evaluate the difference between the first order formalism and the second order formalism 
for our case. We begin with the assumption that the proper time dependence of the shear viscosity 
is given by (|2~9|) . Substituting (|2~9|) to (|B.3j) . we find 

$ = ^|l + (l + /?)^ + (l + /3)(2 + /3)^) 2 + ...| 

+ constant x e~ T ^ Tn . (B.4) 



This means that the second order formalism approaches to the first order formalism when 

- « 1. (B.5) 
r 

Therefore, our analyses based on the first order formalism (|2.8f) with assumption ()2.9j) and the late 
time approximation are self-consistent. The condition (JB.5JI also agrees with our basic assumption 
that the microscopic time scale t v is short enough comparing to the macroscopic time scale so that 
the hydrodynamic description is valid. Small value of r n also matches the fact that our fluid consists 
of strongly interacting particles. 



C Consistency of 77 = 770/V with 77 ~ T 3 



Let us check self-consistency of our assumption H2.13j) . Starting with 77 = r/o/r, the energy density is 
given as 

P(r) = -%(1-2 1 ), (C.l) 
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where 



7 = ^73- (C2) 



POT 

The relationship ij oc T 3 oc p 3 / 4 makes further corrections to the shear viscosity like 

??=-{! + 0(7)}, (C3) 
r 

and the 0(7T _1 ) correction in (|C.3|) makes further corrections to the energy density recursively. 
However, all such corrections are at the higher order of 7 and we can neglect them if 7 <C 1. Let us 
define our approximation precisely: 

• We consider only the region of 7 <C 1. 

• We consider rjT to the order of 1 and /rr 4//3 to the order of 7. In other words, we consider only 
to the order of 770- 

The above makes our framework to be self-consistent. 11 Note that the positive energy condition for 
(|C.1|) is also guaranteed by the above approximation. 
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